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Abstract. This talk reports on recent progress toward the semiglobal study of asymp¬ 
totically flat spacetimes within numerical relativity. The development of a 3D solver 
for asymptotically Minkowski-like hyperboloidal initial data has rendered possible the 
application of Friedrich’s conformal held equations to astrophysically interesting space- 
times. As a hrst application, the whole future of a hyperboloidal set of weak initial data 
has been studied, including future null and timelike inhnity. Using this example we 
sketch the numerical techniques employed and highlight some of the unique capabilities 
of the numerical code. We conclude with implications for future work. 


The modern treatment of gravitating isolated systems allows to study asymptotic 
phenomena by local differential geometry. The principal underlying idea, pioneered 
by Penrose [1], is to work on an unphysical spacetime obtained from the physical 
spacetime by a suitable conformal compactification. Friedrich has extended this 
idea in a series of papers [2] to the level of the field equations by rewriting Einstein’s 
equations in a regular way in terms of equations for geometric quantities on the 
unphysical spacetime Ai. A metric gab on At which is a solution of the conformally 
rescaled equations gives rise to a physical metric cjab = ^~^gab, where the conformal 
factor n is also determined by the equations. The physical spacetime Ai is then 
given hy Ai = {p G Ai \ f2(p) > 0}. These “conformal held equations” render 
possible studies of the global structure of spacetimes, e.g. reading off radiation 
at null inhnity, by solving regular equations. It is natural to utilize everywhere 
spacelike slices S* in Ai which cross null inhnity. On Ai the corresponding slices 
Tit are similar to the hyperboloid — y‘^ — = lA in Minkowski spacetime, 

and are therefore usually referred to as hyperboloidal slices. They are only Cauchy 
surfaces for the future domain of dependence of initial slice of AT, we therefore call 
our studies semiglobal. 

The conformal held equations address a number of problems with the numerical 
treatment of isolated systems in general relativity: Radiation quantities can only 
be dehned consistently at null inhnity (JT’*’). Artihcial outer boundaries cause am- 


biguities and stability problems. Resolving the different length scales of radiating 
sources and the asymptotic falloff is numerically difficult. Our compactified grid is 
allowed to extend beyond the physical part of the unphysical spacetime, the bound¬ 
ary thus can not influence the physics of a simulation. No artihcial cutoff at some 
large distance is required to keep the grid hnite and there is no necessity to treat 
very large length scales (dominating the asymptotic falloff) along with variations on 
small scales. Including in the computational domain enables straightforward 
extraction of radiation quantities involving only well-dehned operations without 
ambiguities. The symmetric hyperbolicity of the implemented formulation of the 
conformal held equations guarantees a well-posed initial value formulation. 

The rest of this article is organized in three parts: First the current technology 
of the solution of the constraints will serve as an example of how some technical 
problems which are particular to the conformal held equations have been solved 
successfully. After outlining the evolution algorithm we discuss the evolution of 
weak data as an example of a situation for which the usage of the conformal held 
equations is ideally suited: the main difficulties of the problem are directly ad¬ 
dressed and solved by using the conformal held equations. These two topics sum 
up the work of Hiibner on 3D numerical relativity with the conformal held equations 
(see [3], [4] and references cited therein). Finally, we discuss future perspectives 
for handling strong held situations, where a number of problems appear which are 
not directly addressed by the conformal held equations, but where we believe that 
their use will prove benehcial. 

The constraints of the conformal held equations (see Eq. (14) of Ref. [5]) are 
regular equations on the whole conformal spacetime Ai. However, they have not 
yet been cast into some standard type of PDE system, such as a system of elliptic 
PDEs. One therefore resorts to a method where one hrst obtains data for the 
Einstein equations - the hrst and second fundamental forms hah and kat induced 
on S by Qat with corresponding Ricci scalar and covariant derivative denoted by 
and Va. After extending this subset of data from S to S the data are then 
completed by using the conformal constraints. Here we restrict ourselves to a 
subclass of hyperboloidal slices where initially kab is pure trace, kab = \habk. The 
momentum constraint V^kab -Vak = 0 then implies k = const ^ 0. We always set 
/c > 0. In order to reduce the Hamiltonian constraint 

‘'^R + P = 

to one elliptic equation of second order, we use the standard Lichnerowicz ansatz 

hab D (j) hab- 

The free “boundary dehning” function D is chosen to vanish on a 2-surface S - 
the boundary of S and initial location of - with non-vanishing gradient on S. 
The topology of S is chosen as spherical for asymptotically Minkowski spacetimes. 
Let hab be a metric on S, with the only restriction that its extrinsic 2-curvature 
induced by hat on S is pure trace, which is required as a smoothness condition [6]. 



With this ansatz hab is singular at S, indicating that S represents an inhnity. The 
Hamiltonian constraint then reduces to the Yamabe equation for the conformal 
factor 0: 


4 V“ Va0 - 4 ) (Va0) 


+ 2QAQ - 3(V“fi)(Vafi) 



This is an “elliptic” equation with a principal part which vanishes at the boundary 
S for a regular solution. This determines the boundary values as 9(V“r2)(VaO) = 
Existence and uniqueness of a positive solution to the Yamabe equation 
and the corresponding existence and uniqueness of regular data for the conformal 
held equations using the approach outlined above have been proven by Andersson, 
Chrusciel and Friedrich [6]. 

The complete set of data for the conformal held equations is obtained from the 
conformal constraints via algebra and diherentiation. This however involves di¬ 
visions by the conformal factor which vanishes at S. In order to 

obtain a smooth error for this operation, the numerically troublesome application 
of I’Hospital’s rule to g = f /Q is replaced by solving an elliptic equation of the 
type — klf) = 0 for g. For the boundary values kl'^g — f2/ = 0, the 

unique solution is g = f /kl. For technical details see Hiibner [3]. The Yamabe 
equation for 0 and the linear elliptic equations arising from the division by are 
solved by pseudo-spectral collocation (PSC) methods (see e. g. [7]). Fast Fourier 
transformations converting between the spectral and grid representations are per¬ 
formed using the FFTW library [8]. Nonlinearity in the Yamabe equation is dealt 
with by a multigrid Newton method (for details and references see [3], the resulting 
linear equations are solved with the AMG library [9], which implements an alge¬ 
braic multigrid technique. The PSC method restricts the choice of gridpoints to 
be consistent with a simple choice of basis functions. Thus S is required to be a 
coordinate isosurface and the elliptic equations are solved in spherical coordinates. 
The polar coordinate singularities are taken care of by only using regular quan¬ 
tities in the computation, i.e. Cartesian tensor components, and by not letting 
any collocation points coincide with coordinate singularities. In order to extend 
the initial data to the Cartesian time evolution grid on the extended hyperboloidal 
slice S, the spectral representations are used, which dehne data even outside of 
the physical domain where the constraints have been solved. The constraints will 
be violated in the unphysical region, but since this region is causally disconnected 
from the physical interior by 77^, the errors in the physical region converge to zero 
with the discretization order. For numerical purposes the coefficient functions of 
the spectral representation are modihed in the unphysical region. 

Time evolution of the conformal held equations (in particular here this is the 
system Eq. (13) of Ref. [5]) is carried out by a 4th order method of lines with stan¬ 
dard 4th order Runge-Kutta time evolution. Spatial derivatives are approximated 
by a symmetric fourth order stencil. To ensure stability, dissipative terms of higher 
order were added consistently with 4th order convergence, as discussed in section 



6.7 of Ref. [10]. In the unphysical region of M near the boundary of the grid, a 
“transition layer” is used to transform the conformal Einstein equations to simple 
advection equations. A trivial copy at the outermost gridpoint yields a simple and 
stable outer boundary condition. 

The main result so far, obtained by Hiibner, is the evolution of weak data which 
evolve into a regular point representing future timelike inhnity. This result illus¬ 
trates a theorem by Friedrich [11], who has shown that for sufficiently weak initial 
data there exists a regular point of M. The complete future of (the physical 
part of) the initial slice was reconstructed in a finite number of computational time 
steps, and the point has been resolved within a single grid cell. For these evolu¬ 
tions very simple choices of the gauge source function have proven sufficient: a zero 
shift vector was used, the lapse was set to = v^duFh, and the scalar curvature of 
the unphysical spacetime was set to zero. The initial conformal metric is chosen in 
Cartesian coordinates as 


/ 1 + (3,2 ^ 2|/2) 0 0 \ 

h= I 0 ^ 

V 0 0 1/ 

and the boundary defining function as = \{f — +1/^ + ^^)) • The extraction 

of physics is largely based on the integration of geodesics concurrently with the 
evolution. This is carried out with the same 4th order Runge-Kutta scheme that is 
used in the method of lines. Null geodesics along are used to construct a Bondi 
system at and to compute the news function and Bondi mass. To illustrate the 
results we show the behavior of geodesics in the numerically generated spacetime. 
First, in Fig. 1 we show three timelike geodesics originating with different initial 
velocities at the same point (xo,yo,zo) = ( 2 ^’ 2 ^’ 2 ^) meeting a generator of 
77 ^ at Fig. 2 shows the oscillations induced by gravitational waves in the zero 
velocity geodesic starting out at (xo,yo,xo) (Fig. (6) of Ref. [4] shows the same 
plot using coordinate time instead of proper time). 

Future efforts will focus on developing the techniques to evolve strong data with¬ 
out symmetries, e.g. to describe generic black hole spacetimes. Here many problems 
are not yet solved, in particular issues associated with choosing the gauge source 
functions, the treatment of the appearance of singularities, and the limitations of 
computer resources for 3D calculations. These well known problems plaguing 3D 
numerical relativity will have to be addressed and solved in the conformal approach 
in order to harvest its benefits for studying the (semi)global structure of spacetimes. 

A first step toward the study of black holes with the conformal field equations is 
to obtain suitable initial data. In the Cauchy approach topologically nontrivial data 
are easy to produce by compactification methods (see e.g. [12], [13] and references 
cited therein) where the topology of the computational grid is not influenced by 
the number of asymptotic regions. In the hyperboloidal case the adding of 77’'“’s as 
suggested e.g. in [3] does change the computational domain and leads to technical 
problems, for which Hiibner has suggested solutions [3]. An alternative is to use 
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FIGURE 1 . Three timelike geodesics (starting at x = meet a generator of (starting 
at X = 1) at future timelike infinity i'^. 
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FIGURE 2 . The differences x — y, x — z, and y — z along a timelike geodesic show the spacetime 
distortion due to gravitational waves. In flat space they would vanish for symmetry reasons. The 
cutoff in time was chosen to better resolve the early time structure. 





regular initial data containing one or more apparent horizons. Such data could be 
produced by parameter studies with the current code and would in some sense be 
more physical, since they do not require the existence of “eternal” black holes, and 
in principle should allow initial data which have apparent horizons and singularities 
in the fnture, bnt not in the past. An important question is whether such data 
are qualitatively different from data with “topological” black holes in the region 
outside of the event horizon. In order to handle the singnlarities inside of black 
holes with the code, a strategy snccessfully employed by Hiibner in the spherically 
symmetric case was to handle floating point exceptions and then continne past the 
formation of singularities as permitted by causality [14]. However, the probably 
chaotic structure of the singularity and the difficulties of evolving with a time step 
as large as permitted by causality may turn ont prohibitive in 3D. Excision of 
the singnlarity or a modification of the eqnations inside the black hole (analogous 
to the way they are changed in the nnphysical region beyond J^') are possible 
alternatives. 
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